On violation of the Robinson's damping criterion and enhanced cooling of ion, 
electron and muon beams in storage rings 
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Limits of applicability of the Robinson's damping criterion []| and the problem of enhanced cooling 
of particle beams in storage rings beyond the criterion are discussed. 
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I. INTRODUCTION 

Losses of energy by particles in storage rings caused 
by friction forces can lead to damping of their betatron 
and synchrotron oscillations and to a decrease of the six- 
dimensional phase space volume occupied by particles 
in a beam (cooling of a beam). Friction can be caused 
by: 1) radiation reaction in static external fields, 2) ra- 
diation reaction in time dependent external fields (back- 
ward Compton and Rayleigh scattering of laser light), 
3) ionization and radiation processes in material targets 
located in the vacuum chamber of the ring. Special inser- 
tion devices can be used to introduce friction forces. The 
physics of damping in the longitudinal plane is different 
from that in the transverse one. 

Usually a particle with a higher than ideal energy loses 
more energy than the ideal particle and a particle with 
lower energy loses less energy. The combined effect is 
that the energy difference between non-ideal and ideal 
particles is reduced. A decrease of the particle energy 
deviation from the ideal leads to a decrease of its longi- 
tudinal phase deviation from the synchronous phase. 

In the transverse plane, the loss of particle energy leads 
to a loss of its longitudinal as well as transverse momen- 
tum, since the particle performs betatron oscillations. 
The total momentum loss is, however, replaced in the 
cavity only in the longitudinal direction. The combined 
effect of the energy loss and the replacement of the en- 
ergy loss in accelerating cavities leads to a net loss of 
transverse momentum. 

The total amount of damping in all degrees of freedom 
is limited by the amount of the energy loss. There is a 
correlation of damping decrements of different planes de- 
termined by Robinson's damping criterion Q-pJ. This 
criterion limits the rate of particle cooling in storage 
rings. Below, we would like to pay attention to the limits 
of applicability of this criterion, its violation in schemes 
of selective interaction of particle beams with targets. 

We start by reviewing Robinson's damping criterion. 



II. LIMITS OF APPLICABILITY OF THE 
ROBINSON'S DAMPING CRITERION 

The motion of particles in storage rings is described by 
their deviations from the ideal orbit in transverse radial 
x, vertical y directions and by deviation ip = ip — ip s of 



the particle phase from the synchronous one in a curvi- 
linear coordinate system (x,y,ip). In a linear approxi- 
mation deviations (x,y,y>) are described by linear second 
order differential equations. We can write the equations 
in the form of a system of six linear first order differen- 
tial equations for a six-dimensional coordinate vector u 
with components (x,x',y,y',ip,A.e), where x' — dx/ds; 
y' = dy/ds; s, the longitudinal coordinate of a particle 
along the ideal (reference) orbit; Ae — e — e s , the devia- 
tion of the particle energy from synchronous one. In the 
matrix presentation: 



du(s) 
ds 



Q(s)u(s) 



(1) 



where Q(s) = is a six-order matrix with com- 

ponents qji(s) — 1...6). The Eq.([Q| has six lin- 
ear independent solutions Uj(s) with components Uji(s) 



(Ujl 



Vj, U j2 



j , Ujs — yj , 



u j6 = Asj). The 



solution of the Eq. has the form u(s) = U(s) ■ u(0), 
where U(s) = \\uij(s)\\ is a transfer matrix; u(0), the 
initial vector. The determinant of this matrix is a Wron- 
skian W(s) = \U(s)\. It represents the six-dimensional 
volume of the polyhedron in the phase space occupied 
by the beam. The values dW(s)/ds = SpQ ■ W(s), 
W(s) = W(0)-exp(J SpQ-ds) ~ W(0)-exp(< SpQ > -a), 

where SpQ = Y^=iQjj> W(0) is the initial Wronskian, 
sign <> denotes averaging. This is the Jacobean for- 
mula, @- On the other hand, Uji(s) ~ exp(ai • s) and 
the rate of change of the 6-dimensional volume of the 
polyhedron ~ exp[2^3 °H s ]-> where on = a x ,a y ,a £ are 
fractional averaged damping decrements. Therefore 



a 6D = 2 oii =< SpQ > 



(2) 



SpQ is determined by the diagonal elements of the ma- 
trix | |<Zji(s)| |. In the transverse plane, the particle mo- 
mentum loss does not lead to a change of the direction of 
the momentum and position of the particle (qn = 933 = 
955 = 0). Acceleration of the particle changes the direc- 
tion of the momentum on the value |Ap|/|p] = P s -s/c-s s . 
It leads to matrix elements (722 = 944 = —Ps/c-e s , where 
P s is the average rate of particle energy loss; e, the par- 
ticle energy; subscript s refers to the reference orbit; c, 
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the velocity of light. The rate of change of the parti- 
cle energy is de/dt = — (dP/de)\ s ■ Ae + (dP r f /dip)\ s ■ <p 
and matrix element ^66 = —{dP/c ■ de)\ s . Substitution 
of diagonal matrix elements to leads to generalized 
Robinson damping criterion 



with rapidly increasing linear dependant spectral inten- 
sity in the frequency range corresponding to the ion en- 
ergy spread er £j o- In this case P/P s ~ (e — s s +^e,o)/<^e,o 
(-a s ,o <£-£ s < a- e ,o), {dT/de)\ s = P s /a e , > P a /e s , 

Ji,max — E-s I^Ue.Q-! 



3 l 



IP 

C 



2<7 f 



, 



The proof of the Robinson's damping criterion was 
reduced to application mathematical Jacobean formula. 
Non-diagonal matrix components responsible for the 
beam dynamics of particles in a lattice was not used. 
Two diagonal components responsible for damping in the 
transverse plane are determined by the average power of 
the particle energy loss and one diagonal component re- 
sponsible for damping in the longitudinal plane is deter- 
mined by the partial derivative of the power energy loss. 
The value cxqd in J5| determine the rate of damping of 
the 6-dimensional phase space volume (emittance) occu- 
pied by the beam (cooling). Coefficients a e and ago can 
be both positive and negative 01- El 

If there is no coupling between radial x and vertical y 
planes in a storage ring, the direct calculations can be 
performed separately for vertical and longitudinal damp- 
ing coefficients: 



.1^ 

2c e s 



1 dP 
2~c ~de ' 



(4) 



The radial decrement follows from Eq. 



(7) 



The value t 6 d = — l/ca^p determine the damping 
time of the six-dimensional phase space volume of the 
beam. In case of radiation of particles in external fields 
the value tqo > (cooling). At the same time in this 
case the damping in one plane and antidamping in an- 
other one is possible . In case of ionization losses of en- 
ergy by muons at small energies the damping time of the 
six-dimensional phase space volume of the beam tqd < 
(heating) 0. 

Next conditions were used to prove the Robinson's 
damping criterion: 1) cooling in the radio frequency (RF) 
bucket, 2) linearity of the system 3) stationary condi- 
tions. Violation of these conditions can lead to the vi- 
olation of the criterion, the concept of decrement (non- 
exponential damping) and fast cooling. 

Below we consider schemes of enhanced emittance ex- 
change and six-dimensional enhanced cooling of particle 
beams beyond conditions used to prove the Robinson's 
damping criterion. 

The term enhanced we refer to the emittance exchange 
and cooling times determined approximately by the equa- 
tion J7J) and less. During this time the energy losses of 
particles of a beam are about the energy spread of the 
beam. 



1 

2c 



£s 



dP 
57 1 



dP 

de 1 



(5) 



Damping times n = — 1/ con are limited by Robinson's 
damping criterion @ by the value 



(6) 



•hP s 



where Ji is determined by the dependence P(e). All 
decrements must be Ji > 0. Every decrement Ji < 
Ji,max — 2 if radiation in external fields is emitted 
(synchrotron, undulator, backward Compton scattering 



one, P 



In case of radiative ion cooling based 



on backward Rayleigh scattering in the homogeneous 
laser beam having uniform spectral distribution the value 
•h < Ji,max = (3 + 2D) /2(1 + D) , where D is the satura- 
tion parameter (P ~ D/(l + D), D ~ e) 0- For ioniza- 
tion muon cooling (dP/de)\ s is rapidly decreasing with 
the energy for e M < 0.3 GeV, but is slightly increasing for 
e M > 0.3 GeV (Ji yma x ~ 1) @- @ The partial derivative 
(dP/de)\ s can be very high only in case of laser cool- 
ing of ion beams by homogeneous broadband laser beam 



III. ENHANCED DAMPING SCHEMES BASED 
ON SELECTIVE INTERACTION OF PARTICLE 
BEAMS WITH MOVING TARGETS 

In this section two schemes of enhanced damping of 
particle beams in the longitudinal and transverse planes 
are considered when external selective interaction of par- 
ticles with moving targets is used ^(j " S3 ■ The targets 
(material, laser beam, undulator) are installed in the vac- 
uum chamber of the storage ring and can be displaced 
in the radial direction. The RF acceleration system is 
switched off and there is no energy dependence of the 
average rate of the particle energy loss. Material tar- 
gets can be used for muon beams and laser targets can 
be used for ion and electron beams (laser beams propa- 
gate in the opposite direction to ion or electron beams) 1 . 



1 Stabilization schemes allow parallel displacement of t he laser 
beam relative to the ion beam to better than 10 /im ll.'J . In 
reality, closed orbits can be moved in the direction of the target 
instead of moving the target. A kick, decreasing of the magnetic 
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The interaction region (IR) in this case changes its ra- 
dial position in the particle beam during the damping 
process. Broadband laser beams have to be used for ion 
beam damping in order for all ions, independent of their 
energy, to interact with laser beams. 

The enhanced damping schemes are presented in Fig- 
ures 1 and 2. In these Figures the axis "z" is the un- 
wrapped reference orbit of the storage ring, L/ and a 
are the length and the width of the targets in the inter- 
action region, respectively. The transverse positions of 
targets are displaced along the axis "x" with the veloci- 
ties i>Ti > and vt 2 < relative to the reference orbit. 
1,2,3,.. are the particle trajectories before and after the 
energy loss in the target. 0\, O2, O3, . . . are the loca- 
tions of the particle closed orbit after 0,1,2,3, . . . energy 
loss events. Parts of one period of the particle trajectory 
at the targets T\ and T2 are shown before and after every 
event of the particle energy loss in the targets. 




FIG. 1: Scheme of enhanced damping of a particle beam in 
transverse plane. The evolution of the amplitude of betatron 
oscillations is shown for the target at rest. 



A target T\ can be used for damping of betatron oscil- 
lation amplitudes of particle beams (see Fig. 1). At the 
initial moment it overlaps a small internal part of the par- 
ticle beam in the radial direction in the straight section of 
a storage ring with non-zero dispersion function. The de- 
gree of overlapping of particle beam and target is changed 
by moving the target position uniformly with some veloc- 
ity from inside in the direction of the particle beam. 
First, particles with largest initial amplitudes of betatron 
oscillations interact with the counter-propagating target. 
Immediately after the interaction and loss of energy, the 
position and direction of momentum of a particle remain 
the same, but the closed orbit is displaced inward in the 
direction of the target |^ . The radial coordinate of the 
closed orbit and the amplitude of betatron oscillations 



are decreased by the same value owing to the disper- 
sion coupling. After every interaction, the position of 
the closed orbit approaches the target more and more, 
and the amplitude of betatron oscillations is reduced. 

The closed orbit moves with a velocity x n , which de- 
pends on the particle's amplitude of betatron oscillation 
and the distance between the orbit and the target, since 
these values determine the probability of collision of a 
particle and a target (see below). It reaches a maximum 
velocity lijjml (x v in < 0) if the depth of penetration of 
the particle closed orbit in the target is greater than the 
amplitude of betatron oscillations of this particle 2 . Af- 
ter this, the amplitude of betatron oscillations will stay 
constant since the particle will interact with the homo- 
geneous target every turn for both positive and negative 
deviation from the orbit with equal probability. When 
the target reaches the closed orbit corresponding to par- 
ticles of maximum energy and the depth of penetration 
of the particle closed orbits in the target is greater than 
the amplitude of betatron oscillations of these particles, 
it must be returned to its previous position. All parti- 
cles of the beam will have small amplitudes of betatron 
oscillations and increased energy spread. 

Particles with high amplitudes of betatron oscillations 
start interacting with the target first. Their interaction 
time is longer and hence the decrease of amplitudes of 
betatron oscillations is greater. 




FIG. 2: Scheme of enhanced damping of a particle beam in 
longitudinal plane. The evolution of the amplitude of beta- 
tron oscillations is shown for the target at rest. 

A target Ti can be used for damping of the energy 
spread of particle beams (see Fig. 2). The radial target 
position must be moved uniformly with a velocity \vt 2 I > 
\&n in I from outside in the direction of the particle beam 
being cooled. At the initial moment, the target overlaps 
only a small part of the particle beam. The degree of 
overlapping is changed in such a way that particles of 



field in bending magnets of the storage ring, a phase displace- 
ment or eddy electric fields can be used for this purpose. 



2 The value x ri 



Sx v /T < 0, T = 1// is the revolution period. 
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maximum energy first come into interaction and then 
particles of lower energy. When the target reaches the 
orbit of particles of minimum energy it must be switched 
off and returned to the previous position. In this case 
the difference in duration of interaction and hence in the 
energy losses of particles having maximum and minimum 
energies will be large. As a result, all particles will be 
gathered at the minimum energy in a short time. 

The evolution of the amplitude of betatron oscillations 
and position of the closed orbit of a particle is determined 
by its energy loss in the target. Contrary to the previous 
case the energy loss of the particle leads to increase of 
the amplitude. But in this case the radial velocity of 
the target displacement \vt 2 \ > \x V in\ and hence after a 
short time the particle closed orbit will be deepened into 
the target to a depth greater than the amplitude of its 
betatron oscillation and the increase will end. 



A. Theory of interaction of particle beams with 
transversely moving targets 

Below we will neglect the energy dependence of the av- 
erage rate of the particle energy loss in the target, the 
emission of SR by particles in the bending magnets of 
a storage ring, assuming the RF system of the ring is 
switched off, targets are homogeneous and have sharp 
edges in the radial directions. We assume the jump of 
closed orbits of particles caused by their energy loss in 
the target is much less than amplitudes of betatron os- 
cillations. 

In a smooth approximation, the motion of a particle 
relative to the position x v of its closed orbit is described 
by the equation xp = Aq cos ipp 7 where xp = x — x v is 
the particle deviation from the orbit; x, its radial coordi- 
nate; ipp = LUp(t~to)+ippo, phase of the particle betatron 
oscillations; (ppo and to, the initial phase and time; Aq 
and uip = 2i:v x j ', the initial amplitude and the frequency 
of betatron oscillations, respectively; u x , the tune of be- 
tatron oscillations; /, the particle revolution frequency 
4J. If the coordinate xpo and transverse radial velocity 
of the particle xpo = — A^ivp sin^o correspond to the 
moment to of change of the particle energy in a target, 
then the amplitude of betatron oscillations of the particle 

before an interaction is Ao = \J x 2 ^ + i^ /u>p. After the 

interaction, the position of the particle closed orbit will 
be changed by a value 5x 71 = r\ x (3 2 (Se/e), where i] x is 
the dispersion function of the storage ring; /3, the rela- 
tive velocity of the particle; 6e, the change of the particle 
energy. The deviation of the particle relative to the new 
orbit will he xpo — 5x Vl and the direction of the particle 
velocity will not be changed. The new amplitude will 

be A\ — ^J{xpo ~ 5x v ) 2 + x^g/uip and the change of the 

square of the amplitude 



In the approximation \6x n \ <C \xpo\ < Aq, the value 



6A = -^6x Tl 



(9) 



The rate of change of the particle betatron oscillation 
amplitude is maximal when \xp§\ = A. 

The velocity of a particle orbit x v depends on the prob- 
ability W of collision of particle and target in one turn. 
This probability is determined by the distance Xt x 2 — x v 
between the edge of the target and the particle closed 
orbit and on the amplitude of betatron oscillations. Par- 
ticles interact with the target at every turn (W = 1) 
and their radial velocity reaches the maximum value x^ i n 
when the orbit enters the target at a depth larger than 
the amplitude of betatron oscillations. In the general 

The deviation of a particle from its closed orbit at the 
location of the target takes on values xp n = Acostpp n , 
where tpp n = 2Tty x n + fpo, n = 1,2,3,... . Different 
values xp n in the region (—A, A) will occur with equal 
probability if v x ~ p/q, q S> 1 (tune is far from forbidden 
resonances), p and q are integers. The collision of targets 
and particle beams occur when xp n < Xt x —x v and xp n > 
xt 2 ~ Xr). These conditions are valid when ipp n are in the 
range of phases 2ipx 1 2 , where (pT x — 7r — arccos£i and 
ipr-2 = arccos^, £1,2 = (a>z\ 2 — x n )/A. Particles cross 
through the target at the range of phases 2<^Ti 2 an d cross 
over it at the range of phases 2-7r— 2ipx 1 2 ■ The probability 
can be presented in the form W = ip^ 2 fit and the value 
x-q = fTi,i ■ x n i n /Tr, where x vin = -i] x f3~ 2 (P/e). 

The behavior of amplitudes of betatron oscillations 
of particles is determined by ©. Particles cross 
the target at different coordinates x po in the range 
of phases 2(fT 1:2 - That is why the average rate of 
change of amplitudes dA/dx n = — xpo/ A, where — xpo 

= -(1/Vt 1 , 2 )J 01 1 ' 2 2 A cos ipp n dippn = ±Asincip Tl , 2 and 
sincifiT li2 = sinipTi 2 /VTi 2 1 the- signs + and — are re- 
lated to the first and second targets. We took into ac- 
count that the limits of integration a^a and 61,2 depend 
on the target: a\ = arccos£i = ir — (px lt a-i = 0, 61 = n, 
b2 = <£T 2 ■ 

Thus, the evolution of amplitudes and closed orbits is 
determined by the system of equations 



dA 

dx„ 



OXrq Xrq i n 



(10) 



From equations (|10fl and the expression dA/dx v 
[dA/dt]/[dx„/dt], it follows: 



dA 



^71 in 

= -—smcp T . 2 



°T) in 



1 - f 2 

1 SI, 2' 



(11) 



5 A 2 = A 2 



A 2 



-2xpo5x v + (SXrjY 



(8) 



Let the initial closed particle orbits be distributed in a 
region ±a v ,o relative to the location of the middle closed 
orbit and the initial amplitudes of particle radial betatron 
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oscillations Aq be distributed in a region a x .o relative to 
their closed orbits. The spread of closed orbits a v o = 

\V,\(3~ 2 Ko/e). 

Suppose that the initial spread of amplitudes of be- 
tatron oscillations of particles o x _q is identical for all 
closed orbits of the beam. The velocities of the closed 
orbits in a target x n i n < 0, the transverse velocities of 
the targets > and vt 2 < 0. Below, we will use 
the relative radial velocities of the target displacement 
ki,2 = VTy where vt 1 2 = d®T\ 3 /dt. In our case, 

in in < 0, fci < 0, k 2 > 1. 

From the definition of 6,2 we have the relation x n — 
xt 12 ~ 6,2^(6,2)- The time derivative is dx n jdt — 
v t 1i2 ~ [A + £i,2(dA/d£i t 2)]d€i t 2/dt. Equating this value 
to the second term in I|1U|) , we obtain the time derivative 



3gl,2 

dt 



^r] in 



1,2 



■6,2(^/0)6,2) 



(12) 



Using this equation we can transform the first value 
in (fTUfl to the form ±sinc<fT 1 2 (6 2) = (dA/d£i 2) 

(0&,2/0t)/ = (TTfci.2 -'<PTi, a )(dA/dZ 1 , 2 )/[A + 

6,2(^/^1,2)] • ^Ti, a or 9mA/96,2 = ±sin(^ Tl2 / 
[7rfci,2 — (VTi 2 i6,2 sin(^Ti, 2 )]- Substitution Z — ip^ 2 ± 
£1,2 sinc^Ti 2 l ea ds to the solution of the last equation 



A = A Q exp 



'6.2 



±sin y>Ti, 2 d£i,2 

7rfei,2 - (<^t 1i2 ±6,2sin^ Tl2 ) 



V 7rfci, 2 - £(6,2) 



7rfci,2 



7rfci 2 - -£(6,2 



(13) 



where the index corresponds to the initial time io, when 
6,0 = -1,6,0 = 1- Substitution A and dA/d£ 1<2 in (12) 
leads to the time dependence 6,2(0 in the form 



t-t = , nA ° , \lp(h,2, 6,2)|, 

1^77 in I 



(14) 



where ^1,2,6,2) - - d^,a^(Ci,a)/A)[»r*x 
(yjTi 2 = ' = 6,2 sin^^ „)] or, according to (|13|) . 



^(^1,2,6,2) 



6,2 



-y / 7rfci ;2 <i6,2 



£1,2,0 1^1,2 - (VT Ija ±6,2sin(^ Tl 2 )| 3 / 2 



= ± 



V /7r l fc l,2|d6,2 



1,2.0 ~ (VT 1>2 ±6,2sin^ Tl 2 )| 3 / 2 ' 



(15) 



Closed orbits of particles having initial amplitudes of 
betatron oscillations j4q, according to (14), penetrate 
into the target to a depth greater than their final ampli- 



tudes of betatron oscillations A r at a moment t 



f 



^0-0(^1,2,6,2,/)/ \x V in\, where 6,/ = 6(*/) = 1, 6,/ 



6(^/) = — 1- During the interval tf — to, the targets 
T\ t 2 will cross a distance |Axbti 2/ 1 = \ v Ti 2 \{tf ~ to) — 
n\ki,2\\ip(ki,2, ^1,2 j)\A . 

The ratio of the final to the initial amplitude of beta- 
tron oscillations, according to (|13fl . is 



A 



«1,2 



fcl,2 



(16) 



If I6.2 1 < 1, the position of the closed orbit, according 
to the definition of 6,2, can be presented in the form 



x v {t) = x Tl 2 + v Tl 2 (t - t Q ) - A[(6, 2 (t)] • 6,2(t). (17) 

At the moment tf the position of the closed orbit of 
a particle, according to (|17|) . will be determined by the 
equation 



x V:0 - *(fci )2 )A , 



(18) 



where #(fci, 2 ) = -6,2,0 + 6,2,/(AfMo) + ^1,2^(^2, 
6,/)- We have used a condition xt x 2a — x rh o + ^06,2,0- 
The moment to depends on the initial conditions for 
particles Aq, x^.q and targets. If the target has a posi- 
tion xt x 20 = at a moment t = 0, then the target will 
contact particles at the moment to = xq/ vt x 2 , where the 
position of the particle xq — 2^,0 + A)6,2,o corresponds 
to the minimal distance of the particle to the target at 
the moment t — 0. If the last closed orbit of the beam 
will be deepened into the target to a depth greater than 
all amplitudes of betatron oscillation of its particles at 
a moment t > tf imax = max{tf}, the orbits will be dis- 



tributed according to the law x v (t) = x 



f~ 



L n,in 



Xr/,0 ~ %n in 

[t- 7rA)V(|fci, 2 |,6,/)/l 

X >: 

(x Vo + A 6,2,o)6,2,//wt 1>2 ] or 



Xrj,f — Xrj^Q " 



Sl,2 



Aa[i& (kip) - ^(^1,2,6,2,/) 



6,2,0 1 



+ , ' ' ] +X v in ■ t (t > tf im ax)- (19) 
Kl,2 

The duration of the cycle of the enhanced damping of 
particle betatron oscillation amplitudes can be presented 
in the form T x ^\ Ao=<J:Lfi = 2a rh0 /\v Tl 2 \ + (tf - t a ) or 



2<7 n ,Q 7T(T Xi o|0(A;i,2,6,2,/)| 



|fel,2| \x v , 



2o- e ,o 
P 



r-R|0(fcl,2,6,2,/) 



(20) 



where R = <T x . /(T nfi . 

Dependencies \ip(ki,2, 6,2,/)| and \t(fci,2, 6,2,/) deter- 
mined by (jT5|l . ((TH|) and the dependence 1^(^2,6)1 on 
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^2 for k 2 = 1.0, k 2 = 1-1 and k 2 = 1-5 are presented in 
Tables 1-5, respectively. The numerical values in these 
tables permit to analyze the evolution of amplitudes of 
betatron oscillations and closed orbits in time. 

Table 1. The dependencies \ij>(ki, and 



|ftx| 


0.00 


0.01 


0.02 


0.04 


0.08 


0.16 


0.32 


0.64 


1.28 


2.56 


oo 


\i>\ 


oo 


2.42 


1.93 


1.54 


1.21 


0.94 


0.70 


0.50 


0.32 


0.19 







1.00 


1.02 


1.02 


1.00 


0.97 


0.90 


0.79 


0.63 


0.45 


0.28 






Table 2. The dependencies \ip(k2, &,f)\ and ^(fe). 



k 2 


1.0 


1.01 


1.02 


1.05 


1.1 


1.2 


1.4 


1.8 


2.6 




oo 


24.38 


13.8 


6.52 


3.71 


2.10 


1.18 


0.65 


0.35 




oo 


66.3 


36.1 


16.0 


8.50 


4.49 


2.34 


1.20 


0.61 



l + lfrl 

\ki\ 



0.5R[9(k 1 )-ir1>(ki,£ 1 ,f) 



1 



(21) 



According to (fTfi)) . <J x j/(J X fi = l/ e — 0.368, if |fci| = 
|fc 1)e | = l/( e 2 - 1) ~ 0.1565. In this case 



|fl«l — \R<1 — '-ox- 



(22) 



The non-exponential damping time of the particle 
beam in the transverse plane, according to l|20|) and con- 
dition k\ = ki t e can be expressed in the form 



7"s LR«1 



^£[e 2 -l] 



12.8a 



£,0 



p p 

In this method of damping the width of the target 



(23) 



Table 3. The dependence |V>(&2, ^2) |fc 2 =i.o- 
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Table 4. The dependence |*0(fe2, ^2) 
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Table 5. The dependence 1*0(^2, £2)|k 2 =i. 5. 
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B. Enhanced transverse damping 

In the method of enhanced damping of betatron os- 
cillation amplitudes of particles in a storage ring by 
a target T\, the transverse betatron beam size is de- 
creased in accordance with <|16[1 to the value <J x j/a x fi — 
Af/A = ■y/\ki\/\(ki \ + 1). The evolution of closed or- 
bits of particles interacting with the target depends on 
their initial amplitudes of betatron oscillations. First of 
all the target interacts with particles having the largest 
initial amplitudes of betatron oscillations and the low- 
est energies. The position of a closed orbit of these 
particles is changed by the law (|17J) up to the moment 
tf. For the time tf — to the closed orbit will be dis- 
placed relative to its initial position (or the position of 
particles having the same initial orbit but zero ampli- 
tude of betatron oscillations) by the value (|18fl . where 

= 1 + + N) - 

According to l|19|) . the final relative dispersive beam 
size and energy spread are 



a > 2(<7,,j +<t x j). (24) 

C. Enhanced longitudinal damping 

In the method of enhanced damping of the energy 
spread of particle beams in storage rings by a target 
T 2 , the transverse betatron beam size is increased in ac- 
cordance with l|16fl to the value cr x j/a Xl o = Af/Ag — 
\Jk 2 j{k 2 — 1). The evolution of closed orbits of particles 
interacting with the target depends on their initial ampli- 
tudes of betatron oscillations. First of all the target inter- 
acts with particles having the largest initial amplitudes 
of betatron oscillations and the highest energies. The po- 
sition of the closed orbit of these particles x Vl is changed 
by law l|17(l up to the time t = tf. For the time tf — to the 
closed orbit will be displaced relative to its initial position 
(or the position of particles having the same initial orbit 
but zero amplitude of betatron oscillations) by the value 
(HU, where tf(fc 2 ) = Trfc^fca, &,/) " \/k 2 /{k 2 - 1) - 1. 

According to l|19|) . the final relative dispersive beam 
size and energy spread are 

+ 0.5i*[*(fc 2 ) - 7rV>(fc 2 , 6,/) + (25) 

k 2 k 2 

According to (ftT)|) . (T x j/<T x ,a = \fe ~ 1.65, if k 2 — 
k 2 , e — e/(e — 1) — 1-58. In this case 

^|k<i - =±0.37. (26) 

The non-exponential damping time of the particle 
beam in the longitudinal plane, according to (121 )li and 
condition k 2 = k 2>e can be expressed in the form 
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T e|fl«l 



2cr £:0 e-l 1.27(T E) o 



p 



p 



(27) 



In this method of damping the width of the target is 
determined by l|24l) as well. 

According to (jTBjl . JT5J, a rectangle on the plane 
(x^ , A) is transformed to a parallelogram of the same 
area. The four-dimensional phase space volume occupied 
by the beam in this case is not changed in the geometri- 
cal variables {x,x ,y,y , s, x v ). In the canonical variables 
(x,p x ,y,p y , s,p s ) the Liouville's theorem does not work 
and we have a non-enhanced cooling only 3 . 

Friction and the external selectivity of interaction of 
homogeneous moving targets with particle beams in stor- 
age rings do not lead to enhanced cooling if the energy 
independent power loss in the target and the approxi- 
mation \5x n \ <C a x fi are used 4 . In particular, with the 
framework of the above approximation particles are deep- 
ened in the target T2 to the depth larger then their am- 
plitudes of betatron oscillations for many turns, inter- 
act with the target for this time at deviations from the 
closed orbit xp of one sign and, according to re- 
ceive the unwanted increase of betatron amplitudes 5 . If 
particles are deepened in the target to the depth larger 
then their amplitudes of betatron oscillations for one turn 
> Cx,o); they interact with the target at devia- 
tions from their closed orbits of alternate signs. How- 
ever this and similar interactions do not lead to cooling. 
They lead to change the location of the being interacted 
particles from one place to another in the geometrical 
six-dimensional phase space without change of the occu- 
pied by these particles volumes and without overlapping 
these volumes with ones occupied by another particles in 
this space 6 . An ordinary week cooling is observed in the 
phase space determined in canonical variables. Another 
more complicated schemes of external selectivity based 
on moving targets must be used for enhanced cooling. 

The enhanced longitudinal method of laser cooling 
based on internal selectivity in combination with the en- 



3 If we return particle beam to the initial energy state by eddy 
accelerating fields, the transverse phase space volume in the ge- 
ometrical variables is decreased. 

4 Terms "enhanced cooling in one plane" and "enhanced heating 
in another one" used in our previous papers are misleading as 
they concern to enhanced emittance exchange. The term "cool- 
ing" is applied to damping of six-dimensional phase space volume 
determined in canonical variables. 

5 This interaction is similar to interaction of ions and monochro- 
matic laser beams with scanning frequency (see above). However, 
in the last case laser beams overlap being cooled ion beams, in- 
teraction occur at the deviations x@ from their closed orbits of 
different signs, do not lead to increase of amplitudes of betatron 
oscillations or leads to a weak stochastic excitation of betatron 
oscillations if dispersion function is not equal zero. 

6 Another particles previously located in the being occupied vol- 
ume can not stay this volume as they are forced to interact with 
the same target the same moment. 



hanced transverse method of damping can be effectively 
used for three-dimensional cooling of ion beams (see sec- 
tion IV). An enhanced cooling method based on more 
complicated scheme of the external selectivity and mov- 
ing screens located on the way of laser targets can be 
effectively used for three-dimensional enhanced optical 
cooling of proton, fully stripped ion and other particle 
beams (see section V). 



IV. ENHANCED COOLING OF ION BEAMS 
BEYOND THE ROBINSON'S DAMPING 
CRITERION. 

Below enhanced laser cooling of ion beams beyond the 
Robinson's damping criterion is discussed. Internal ion 
selectivity in the process of Rayleigh scattering of laser 
photons is used. Three examples are considered. 

1) Monochromatic laser beam target with scanning fre- 
quency is used when the RF system of the storage ring 
is switched off (RF buckets, linear dependence -P(Ae) 
and stationary conditions are absent) The laser 
beam overlaps the ion beam. Not fully stripped (elec- 
tronic transitions) or naked (nuclear transitions) ions in- 
teract with the homogeneous counter-propagating laser 
beam at resonance energy, decrease their energy in the 
process of the laser frequency scanning until all their en- 
ergies reach the minimum energy of ions in the beam. At 
this frequency the laser beam is switched off. 

The higher the energy of ions, the earlier they start 
interacting with the laser beam, the longer the time of 
interaction. Ions of minimum energy do not interact with 
the laser beam at all. At the same time the amplitude 
of betatron oscillations is not changed, if the dispersion 
function of the storage ring at the IR is zero. As a re- 
sult, the energy spread is decreased by nonexponential 
low to very small value determined by the width of the 
ion spectral line, the bandwidth of the laser beam and 
the average energy of the emitted photons 7 . The radial 
amplitudes and the emittance will not be changed. The 
damping time is determined by Q. 

2) Ion and broadband laser beams interact in a straight 
section of a storage ring. The laser beam is homoge- 
neous in limits of the ion-laser beam IR and has sharp 
frequency edges. The frequency band of the laser beam 
is sufficient for all ions to interact with the laser beam. 
The RF system of the storage ring is switched off (viola- 
tion of conditions 1,2). The minimum initial energy of 
ions corresponds to interaction with laser beam photons 
of high-frequency edge. The dispersion function of the 
storage ring at the IR is zero. 

In this case ions decrease their energy until all their 



7 Exponential damping leads to decrease of the beam dimension 
e ~ 2.7 times for one damping time while non-exponential damp- 
ing leads to much greater decrease and much faster cooling. 
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energies reach the minimum energy of ions in the beam, 
the radial amplitudes and emittance will not be changed. 

3) Ion and broadband laser beams interact in a straight 
section of a storage ring. The laser beam is homoge- 
neous in limits of the ion beam and has sharp frequency 
edges. The RF system of the storage ring is switched 
on. The synchronous energy of ions corresponds to inter- 
action with photons of high-frequency edge of the laser 
beam. The spectral intensity of the laser beam is linearly 
decreased from a maximum at the low-frequency edge to 
zero at the high-frequency edge. 

In this case a discontinuity in the rate of energy loss 
was introduced: ions with an energy more than the 
synchronous energy interact with the laser beam and 
ions with less energy do not. A synchronous ion does 
not lose energy. There is no friction and antidamp- 
ing of synchrotron oscillations at the energy e < e s 
and there is damping at the energy e > e s (violation 
of conditions 1, 2). The power of the scattered radi- 
ation depends on the ion energy according to the law: 

P = Pynax [(e - E s )/(T Ej0 \ at £ s < £ < £ s + cr £ ,o and P = 

at e < e s , s > s s +er e .o. The minimum damping time will 
be determined by if we accept in J7| P s = P-max- The 
radial amplitudes and emittance will not be changed. 

Damping time for enhanced ion cooling J7J) is 
e/4cr £j o > 10 2 times shorter than damping time for ra- 
diative ion cooling © if scattered powers are the same. 

The considered enhanced methods of laser cooling in 
longitudinal plane are based on internal resonance selec- 
tivity installation-specific for ions. In combination with 
the enhanced transverse damping considered in previous 
section they can be used effectively for three-dimensional 
cooling of ion beams (see Appendix 1) The emit- 

tance exchange through a synchro-betatron resonance 
|22| or dispersion coupling by a wedge-shaped laser tar- 
get can be used for three-dimensional cooling as well by 
analogy with the idea of muon cooling [^|. |23tl24ll25| . 



V. ENHANCED OPTICAL COOLING OF 
PARTICLE BEAMS BEYOND THE ROBINSON'S 
DAMPING CRITERION 

Below a method of enhanced optical cooling of particle 
beams based on external selectivity is considered. In this 
method two identical undulators are installed in different 
straight sections of a storage ring with high-dispersion 
and low-beta functions at a distance determined by a 
phase advance 2pn + 7r for the lattice segment, where 
p = 1,2,3... is a whole number. Undulator Radiation 
(UR) emitted by a particle in the first undulator pass 
through an optical system with movable screens located 
in the image plane of the particle beam. Then this radi- 
ation is amplified and pass through the second undulator 
together with the particle. Screens in the optical sys- 
tem open first the way for UR emitted by particles with 
higher energies and higher positive deviations from their 
closed orbits xp > 0. The beam of amplified UR in this 



case is equivalent to moving prototype of the target Ti 
considered above if definite phase conditions are fulfilled 
in the optical system to inject particles in the second un- 
dulator at decelerating phases. In this case energy losses 
are accompanied by a decrease both energy spread and 
amplitudes of betatron oscillations of particles that is by 
enhanced cooling 8 . After the screen will open images of 
all particles of the beam the system must be closed. Then 
the cooling process can be repeated. Laser and optical 
systems in the considered scheme of cooling are similar 
to optical systems in the scheme of the optical stochastic 
cooling |2t|. 

Another schemes of optical cooling using external se- 
lectivity can be suggested. For example, similar scheme 
can use first undulator and even number of undulators in- 
stalled in straight sections of a storage ring at distances 
determined by a phase advance (2p+ l)ir between neigh- 
boring undulators. In this scheme deviations of particles 
in undulators "i" and are xp i — —xp i+1 and that 

is why the decrease of the energy of particles, according 
to (J9J, do not lead to change of their betatron ampli- 
tudes and leads to cooling of the particle beam in the 
longitudinal plane. 

The wavelets of UR emitted by a particle in the first 
undulator and amplified in the optical amplifier inter- 
act efficiently with the particle in the second undulator 
and do not disturb trajectories of another particles if the 
average distance between particles in longitudinal direc- 
tion is less than the length of the UR wavelets KX, where 
K is the number of the undulator periods; A, the wave- 
length of the UR 9 . The efficiency of cooling is higher, if 
the transverse laser beam dimensions of the wavelets in 
the second undulator are less then the transverse disper- 
sion and betatron dimensions of the being cooled particle 
beam. 

Enhanced optical cooling in a RF bucket can be pro- 
duced as well. The damping time of the order of J7J) can 
be received. 

The considered schemes is of great interest for cooling 
of proton, muon and fully stripped ion beams 10 . 



8 If the phase advance for the lattice segment is 2piz + tt and the 
deviation of the particle in the first undulator x@ > 0, the de- 
viation of the particle in the second undulator xp < 0. In this 
case the amplitude of betatron oscillations of the particle in the 
second undulator is decreased. 

9 The amplified UR do not disturb the next particles (if over- 
lapping occur) in the first approximation and leads to a week 
increasing of their amplitudes of oscillation in the second one 
because of the stochasticity of the phase of the UR wavelet for 
another particles. 

Laser cooling b ased on nuclear transitions has problems with 
low-lying levels 1271 . Optical cooling of heavy ions, on the level 
with optical stochastic cooling, is the most efficient |2S|. In this 
case the emitted power ~ Z 2 , where Z is the atomic number. 
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Appendix 1 



Machine parameters 



Below we consider the next scheme of the enhanced 
laser cooling of ion beams. Moving broadband laser beam 
target T\ with sharp high-frequency edge produces damp- 
ing in the transverse plane. Then the target is stopped 
and produces enhanced cooling in the longitudinal plane. 
The RF system of the storage ring is switched off. Non- 
zero dispersion function at the IR is used. 

Ions having maximal amplitudes of betatron oscilla- 
tions, according to lj2"T]l. will be deepened in the tar- 
get after they loose the energy Se/e = a X} o[^(ki ie ) — 
Trip(ki. e ^ij) + l/ki, e ]/r] x necessary for transverse damp- 
ing of theirs amplitudes. That is why the frequency edge 
of the laser beam must correspond to the minimum ini- 
tial energy of ions decreased on the value 5e and the 
frequency band of the laser target must be 



^L L> ^±Se = 2a^ F{RMel 



where F(R, fc x e ) = 1 + 0.5i?[*(fci e ) - 7n/>0i e , £i /) + 

l/h, e }. 

In this case the transverse damping time is determined 
by (|23JI . The enhanced cooling in the longitudinal plane 
is determined by the damping time 



2<r e .o + Se 2(7 



E,0 



F(R, k 1>e ). 



(29) 



For damping time (|29J) ions will be gathered at the min- 
imum energy corresponding to the high-frequency edge 
of the laser beam. The beam will be cooled both in the 
transverse and longitudinal planes. 

Example 

Below we consider an example for cooling of N-like 
Xenon ions (I 29 'Xe 47+ ) in the RHIC storage ring when 
the transition between the ground state (2s 2 2p 3 ) 4 5 3 / 2 
and the excited state (2s 2 2p 4 ) 4 P3/2 is used. The RF ac- 
celeration system of the RHIC is switched off. Enhanced 
ion beam cooling transversely and fast cooling longitudi- 
nally by a homogenous broadband laser with sharp high- 
frequency edge is used. 



Ion properties 



Degeneracy factors 
Transition energy 
Wavelength 
Oscillator strength 
Natural line width 
Decay length 



9i = 92 = 4 
hw tr = 608.44 eV 
X tr = 2.04- 1CT 7 cm 
fa = 8.9 x 1(T 2 
Auj tr /uj tr = 1.58 • 1CT 6 
2gifi. 2 r e uJt r / 92 = 2 cm 



Circumference 
Relativistic factor 
Energy 

Energy spread 
Relative energy spread 
Emittance 
/?— function 

Dispersion function at the IR 
RMS dispersion beam size 
RMS betatron beam size 



C = 3.8 km 
7 = 97 
e = 1.18 ■ 10 13 eV 
2cr £i0 = 4.72 x 10 9 eV 
2cr £i0 /e = 4 x 1CT 4 
e x = 4 x 10~ 7 m-rad 
X = 1 m 
Vx — 5 m 
2a ri .Q = 2 mm 
a x o = 0.632 mm 



Parameters of the laser beam T\ 



Wavelength 
Intensity 

Saturation parameter 
Bandwidth 
Intraresonator power 
RMS radial beam size 
Rayleigh length 
Interaction length 



\ Ll = 3954 A 
I Ll = 1.16 x 10 6 W/cm 2 
D = 10~ 2 
Alu Li /lu Li =9.7-10- 4 
P Ll = 2na 2 Li I Ll ~ 0.45 MW 
o~l 1 — 2.5 mm 
z R = ^a\J\ Ll ~ 20.8 m 
hnt — 10 m 



The required power of the laser beam should be feasi- 
ble to obtain with a free electron laser in an intracavity 
configuration Q 11 or with a conventional laser using a 
high-finesse (~ 10 4 ) optical resonator |3ll] 12 . 

In this case the saturation intensity I sa t — [.91/(51 + 
g 2 )]{2^chio tr / 1 2 \ z tr ) (Alu l /lu l ) ~ 2.87 x 10 8 W/cm 2 ; 
the power of the radiation scattered by ion P = 
2i 2 l mt P L o R /(l + D)CSeff = 8.79 • 10- 10 W = 5.49 • 10 9 
eV/sec ~ 6.95 • 10 4 eV per turn; the effective area of the 
interaction region of laser and ion beams S e f t = 2n{o L + 
ct 2 0) = 0.42 cm 2 ; the Rayleigh scattering cross-section 
a R = Trfi2r e \ tr LUL/Auj L = 6.434 • 10 18 cm 2 ; r e = e 2 /mc 2 
0, lifex I = 0.1565, R = 0.632, F(R,k he ) = 2.36. The 
bandwidth of the laser beam is in accordance with 1281) . 

For these parameters of ion and laser beams, the ve- 
locity of the closed orbit of an ion in the laser beam 
in in = VxP/s — 0.233 cm/sec, the average energy of 
the scattered photons is < Two s >= hcj tr J = 59.0 keV, 
the velocity of the laser beam vt x — 3.64 ■ 10~ 2 cm/sec, 
a xJ = 0.232 mm, 2a eJ = e{Au L /u; L ) = 1.14 • 10 10 eV, 
2cr„ f — 4.85 mm, 2a f < a = 2&i Jl . 

In this case, according to (|23l) , the transverse damping 
time t x = 5.5 sec. and, according to (|29|l the longitudi- 
nal damping time r 6 ~ 2.03 sec. The damping time of 
the same beam in the method of radiative ion cooling, 
according to ©, is t x — r y ~ t £ ~ 4.3 
~ 10 3 times higher 13 . 



10 3 sec, i.e., 



11 At present such lasers are in operation |21| . |3(1 . 

12 Finesse of 1.9 ■ 10 s was reported near A = 850 nm in |32| . 

13 In this scheme a broadband laser beam overlaps an ion beam, all 
ions interact with the laser beam independently of their energy 
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The degree of the reduction in the dispersion beam 
radius in the non-exponentional longitudinal cooling of 
the ion beam can be very small (cr^. //o^.o !)• It 
depends on the sharpness of the high frequency edge of 
the laser beam and quantum processes 14 . 

Cooling of the beam in the longitudinal plane can be 
performed by a monochromatic laser beam with a scan- 
ning frequency for the time determined by (|29|) , which is 
t s <C t x . It can be located in another straight section 
of the storage ring. Unfortunately, using of the optical 
resonator in this case (laser beam with a scanning fre- 
quency) has problems. 

To increase the degree of damping of betatron oscilla- 
tions in the transverse plane we can use a target moving 
back and forth between outside orbit of the beam to in- 
side one. A phase displacement mechanism or eddy elec- 
tric fields can be used for the recovery of the ion beam 
energy and production of the next cooling cycle. 

Schemes of six-dimensional cooling in the RF bucket 
can be suggested as well. 

In the considered method of ion cooling, the laser beam 
does not have a sharp edge. Since the intensity of the 
beam is changed with the displacement of the ion orbits 
to the laser beam center, the laser beam will be simi- 
lar to the wedge-shaped target. Such targets at rest de- 
crease the amplitudes of ion betatron oscillations in case 
of the target T\ and increase them in case of the target 
T2 in lesser degree than the targets with sharp edges. 
This means that the requirements for sharpness of the 
laser target T% are not strong. A system of lasers can be 
used to increase the sharpness of the effective composite 
laser beam. For example, a third laser beam of the same 
maximum intensity as the first one with the dimension 
<7l 3 ~ 0Li/4 intersecting the ion beam at a third straight 
section of the storage ring with the displacement 2ct_l x rel- 
ative to the radial position of the first laser can be used 
in the example. 

Cooled not fully stripped ion beams can be used in 
New Generation Light Sources |36|. 



beams (proton, ion, muon) based on external selectivity 
is suggested and developed. 

Author thanks A.M.Sessler and Robert Palmer for use- 
ful discussion and acknowledges the support of this work 
by the RFBR under Grant No 02-02-16209. 



VI. CONCLUSION 



We investigated the limits of applicability of the Robin- 
son's damping criterion to the problem of cooling of par- 
ticle beams in storage rings. Theory of the emittance 
exchange based on external selectivity caused by mov- 
ing material targets is developed. New schemes of six- 
dimensional enhanced cooling of ion beams based on in- 
ternal selectivity of ions beyond the criterion are consid- 
ered. A scheme of enhanced optical cooling of particle 



and amplitude of betatron oscillations 0, l33l - l35l . The physics 
of radiative ion cooling is similar to SR damping. 
The energy spread of the being cooled ion beam is limited by 
quantum processes of Rayleigh scattering to the value cr e ,q — 
< hw s >= 59 keV. 
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